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THE IMPLICATIONS OF GUNN-PETERSON TROUGHS IN THE HEII LYMAN-a FOREST 
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ABSTRACT 

Many experts believe that the z ~ 3 Hell Lya forest will suffer from the same saturation issues as 
the z ~ 6 HI Lya forest and, therefore, will not be a sensitive probe of Hell reionization. However, 
there are several factors that make Hell Lya absorption more sensitive than HI Lya. We show that 
observations of Hell Lya and Ly/3 Gunn-Peterson troughs can provide a relatively model-independent 
constraint on the volume-averaged Hell fraction of XHeii.v ^ 0.1. This bound derives from first 
using the most underdense regions in the Hell forest to constrain the local Hell fraction and, then, 
assuming photoionization equilibrium with the maximum allowed photoionization rate to calculate 
the ionization state of nearby gas. It is possible to evade this constraint by a factor of ~ 2, but only 
if the Hell were reionized recently. We argue that Hell Lya Gunn-Peterson troughs observed in the 
spectra of Q0302-003 and HE2347-4342 signify the presence of > 10 comoving Mpc patches in which 
2:hoII,v > 0.03. This is a factor of 20 improvement over previous constraints from these spectra and 
100 times stronger than the tightest constraint on the HI volume-filling fraction from the z > 6 HI 
Lyman forest. 

Subject headings: cosmology: diffuse radiation - intergalactic medium - quasars: absorption lines 



1. INTRODUCTION 

Most experts do not believe that effective optical 
depths of Teff ~ 5 measured in the z « 6 HI Lya 
forest imply that HI reionization was occurring (e.g. 
iBecker et alll2007[ ). Observations of the z « 3 Hell Lya 
forest also find t^s 5 in select regions, and again the 
view of experts is that these observations do not neces- 
sarily indicate tha t Hell reionization is happen ing (e.g. 
iFardal et all [19981 : iMiralda-Escude et al.l [2OO0I) . How- 
ever, there are several key differences that make Hell 
Lya absorption at z w 3 a more sensitive probe of dif- 
fuse Hell than HI Lya absorption at z > 6 is of diffuse 
HI: 1) the abundance of helium is « 14 times smaller 
than that of hydrogen, 2) a photon redshifts through the 
Hell resonance 4 times faster than the HI resonance, 3) 
the intergalactic medium (IGM) is less dense and has 
more voids at z « 3 than z « 6, and 4) the HI Lya 
forest absorption is available for each Hell sightline and 
reveals the sightline's density structure. 

Previous theoretical studies of the Hell Lya for- 
est modeled this absorption as a set of dis- 
crete abs orbing clouds a n d/or focused on th e T^a 
statistic (iGiroux fc ShTilil 119971 : iFardal et all 119981 : 
IMiralda-E scude et al.l l200(J) . This Letter employs better- 
suited techniques to study Hell Lya absorption near sat- 
uration. The discrete-cloud approximation works best in 
overdense regions. The low-density IGM is better ap- 
proximated as a less clumpy continuum of gas that has 
not decoupled from the Hubble flow, and it is the ab- 
sorption in the least dense regions that are the most 
constraining with regard to Hell reionization. Further- 
more, the local measure for the Hell fraction studied 
here is more powerful than global measures like Teff(z), 
for which the effect of Hell reionization is highly model- 
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dependent. Our measure uses HI absorption to tar- 
get the Hell absorption in the most evacuated voids. 
The advantages of using the HI absorption were real- 
ized in many observational studies of the Hell forest 
(iDayidsen et al .11 19961: i Hogan et al.lll997l:rAnderson et all 
ll999HHeap et a l. 2000; Smette et al.ll2002[ ). and we argue 
here that this information can be used to obtain tighter 
constraints on the amount of intergalactic Hell than pre- 
vious analyses have found. 

This study is timely because the Hubble Space Tele- 
scope (HST) reservicing mission installed the Cosmic 
Origins Spectrograph (COS) in May 2009. COS wiU en- 
able absorption measurements at 2.76 < z < zqso for 
Hell Lya and at 3.46 < z < zqso for Hell Ly/3, where 
•zqso is the redshift of the quasar. It is able to achieve 
higher signal-to-noise ratios and higher spectral resolu- 
tions than previous instruments on ~ 4 existing Hell 
Lya sightlines, in addition to providing > 10 new sight- 
lines. Many candidate sightlines hav e been identified as 
promising for COS spectroscopy (Sv phers et al.l[2009l ). 

In this Letter, we assume a flat ACDM cosmology 
with h = 0.71, = 0.046, VL^ = 0.27, ag = 0.8, 
and iHe = 0-24, con sistent with recent measurements 
(|Komatsu et al.ll2009f) . 

2. THE HEII LYMAN-SERIES FOREST 

In contrast to HI Lya Gunn-Peterson absorption at 
z « 6, which saturates for neutral fractions of ~ 10^^, 
the Gunn-Peterson optical depths for a photon to redshift 
through the Hell Lya, Ly/3 and Ly7 resonances are 

0.01 ) \ A J \0.1 
(1) 

where A;, is the gas density in units of the cosmic mean 
and xhcII is the local fraction of helium in Hell.^ Lower 

^ Note that damping wing absorption from Hell in dense systems 
or from intergalactic Hell is not significant. 



^Hcii(",/?,7) = (3.4,0.7,0.2) 



2 




2.8 2.9 3 3.1 3.2 

redshift 



Fig. 1. — HI hya forest spectrum (dotted curve) and Hell Lyo 
forest spectrum (solid curve) of Q0302-003 taken with FORS2 on 
the V LT and STIS on the HST, respectively (Worscck & Wisotzki 
[2001). The dot-dashed curve is the la uncertainty on the STIS 
measurement. Modulo uncertainty in the continuum level, the 
long-dashed curve signifies 98% transmission and the short-dashed 
T = 4 or Si 2% transmission. The cyan highlighted regions de- 
lineate underdense regions (as indicated by the HI data) in which 
the Hell Lyo absorption is saturated, and the dark blue region is 
contaminated by geocoronal hya emission. 

redshift HcII Lya absorption obscures the Hell Ly/3 ab- 
sorption (and both Lya and Ly/3 obscure the Ly7). If 
the foreground Lya absorption is from z < 2.8 (where 
TeS, hya ^ 1), this obscuration should not significantly 
weaken bounds that ignore it. 

The focus of this Letter is on absorption in underdense 
regions and what this absorption implies about the ion- 
ization state of denser, neighboring regions. Approxi- 
mately 13% of the volume has Af, < 0.2 at z = 3, 5% 
has Ah < 0.15 (jMiralda-Escude et al.|[2000[ ). and a larger 
fraction in redshift space. Equation ([T]), combined with 
the knowledge that the least dense regions in the z ~ 3 
IGM have Ab ^ 0.1, suggest that the Hell Lya forest is 
sensitive to HeH fractions in such underdensities of the 
order 1% (and higher series hues to ~ 10%).^ 

In order to achieve such constraints, the locations of 
underdense regions must be known. The HI Lya forest 
reveals this information. The HI Lya Gunn-Peterson 
optical depth is 

Thi (a) ^0.7 a, 1^20, 000 j ^ [ 4 J 

(2) 

where we have assumed photoionization equilibrium with 

* Underdense regions will be expanding faster than the Hubble 
flow such that Equation [T] will be an overestimate (by we estimate 
from cosmological simulations < 50% for At, < 0.15). Accounting 
for peculiar velocities does not change our constraints on xjjcn 
because peculiar velocities have the same effect on both Tjjon a^nd 
TBI- 



Ti2 - the HI photoionization rate in units of 10 s ^ 
-, and a power-law temperature-density relation given 
by T = TqAJ-^ where 7 < 1.6 (IGnedin fc Huilfltol) . 

The term A^''^'"'^ in equation [2] arises from the tem- 
perature dependence of the recombination rate. We 
use the values 7 = 1.3, To = 18, 000 K, and = 
0.8, consistent with observatio ns ([McDonald et al.l[2001l : 
[Faucher-Giguere et al.[ [2008a[ ). is expected to be 

nearly spatially invariant owing to the long mean free 
path (mfp) of hydrogen-ionizing photons and the large 
number of sources within a mfp. 

The smaller the value of Af, that can be reliably located 
in the HI Lya forest, the better one can constrain the 
HcII fraction from the Hell Lyman forest. However, con- 
tinuum fitting tends to artificially remove flux such that 
the flux in the lowest density pixels is set to zero, thereby 
preve nting one from distingui s hing b etween small values 
of Af,. [Faucher-Giguere et al.l ([2008[ ) estimated that con- 
tinuum fitting on average removes 3% of the transmission 
at z = 3. Therefore, a value of A;, — 0.15, which yields 
~ 0-03 or 3% absorption, is approximately the min- 
imum density contrast over that can be discriminated 
at z ^ 3. While a more rigorous derivation of the mini- 
mum Ah would be worthwhile, we employ A^ — 0.15 for 
this study. 

Figure [T] shows the HI Lya spectrum (dotted curve, 
R = A/AA = 1000) and Hell Lya spectrum (thick 
solid curve, R = 800) of the z = 3.29 quasar Q0302- 
003 ([Worseck fc Wisotzkill2006l) . one of the most-studied 
Hell Lya sightlines. Note that for 3.1 < z < 3.2 - a 
region spanning « 100 comoving Mpc (cMpc) ~ t < 4 
(except perhaps near z = 3.1) such that there is no sig- 
nificant Hell Lya transmission. The short-dashed hor- 
izontal line in Figure [T] represents the transmission for 
r = 4. 

Many of the elements with t < 4 in Hell Lya absorp- 
tion admit ~ 100% HI Lya transmission (see the cyan 
highlighted regions in Figure [T]) . These regions should 
have Af, < 0.15 from the above discussion. Equation ^ 
and the limit r < 4 implies that the underdense elements 
with Af, < 0.15 have a;Hcii > 0.008.^ Throughout this 
Letter, we adopt this bound on Xfjoii for the most un- 
derdense, saturated elements in the Hell Lya spectrum 
of Q0302-003, although a more detailed analysis could 
improve upon it. Another well-studied Hell forest sight- 
Hne, the z = 2.89 quasar HE2347-4342 (not shown), has 
a Gunn-Peterson trough at z « 2.85 with no detected 
Hell Lya transmission and provides a similar bound. 

3. HEH PHOTOIONIZATION 

Sometime before z « 2.8, the second electron of inter- 
galactic helium was reionized by a source of ultraviolet 
photons (most probably quasars), and afterward the he- 
lium was kept doubly ionized by the metagalactic radi- 
ation background. If a gas parcel is exposed to a radi- 
ation background with Hell photoionization rate FhoII, 
the Hell fraction as a function of elapsed time At is 

XHcuit) ~ XHoII.oq + (a^HcII.O - 2:HoII,oq) CXp (-At/toq) , 

(3) 

^ We have checked that this bound holds if we instead use the 
R = 40, 000 Keck HIRES HI Lya forest spectrum of Q0302-003. 
Since underdense regions tend to have several hundred km s^^ 
widths, their character does not change with higher resolution data. 



3 




-16.5 -16.0 -15.5 -15.0 -14.5 
LogiorHeii 

Fig. 2. — Hell fraction of a patch in the z = 3 IGM as a func- 
tion of the value of rHoii, assuming photoionization equilibrium. 
The thin and thick solid curves are the Hell fraction at A;, = 0.15 
and Aj, = 1. The short- and long-dashed curves are respectively 
the volume- and mass- weighted HeH fractions, including only ele- 
ments with A[, < 10. The intersection of the solid vertical line (at 
FHeii ~ 6 X 10~^^ with these curves sets the lower bound on 
the respective Hell fraction if regions with A(, = 0.15 have r > 4 
(our constraint) . The dashed vertical line is a plausible future con- 
straint that is 5 times stronger. 

where .thoIi.o is the Hell fraction at At — 0, t^q = 
(rHoii + a71e)"\ 

_ a (T) fie , . 

rHcii + a (T) Tie 

Uf, is the electron density, and ol{T) is the recombination 
coefhcient.^ The z ~ 3 recombination time [a(T)ne]~^ 
is roughly half the Hubble time H{z)~^ at mean density. 
Figure[2]plots XHeii.eq as a function of FHeii for Af, = 0.15 
and Ab = 1. 

4. A NEW METHOD TO ESTIMATE Xheh 

We have shown that a measurement with percent-level 
precision of the HeH Lya forest can be used to set a 
percent-level constraint on the HeH fraction in under- 
dense regions that are opaque to HeH Lya photons. In 
this section, we argue that such a measurement can place 
an even stronger limit on the local Hell fraction in a 
> 10 cMpc region surrounding such underdensities. 

Our argument requires two propositions to hold: 

1. The ionization state of gas elements in a volume 
that is exposed to the same FhcIiIO scales as 
xhoII fx the scaling expected for pho- 
toionization equilibrium.'' 

2. A > lOcMpc region surrounding an underdensity 
in the forest is exposed to approximately the same 
rHcii(i) as the underdensity. 

® Note that A;, appears implicitly in a (T) ne and a is the case 
A recombination coefficient. 

This proportionality omits the A;, factors that appear in the 
denominator of Equation|4]for notational convenience (see footnote 
5). Our calculations include this additional scaling, which increases 
in importance with density. 
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Fig. 3. — Illustration of possible scenarios in the A;, versus 
2;hgII plane at 2 = 3. The solid curves assume equilibrium with 
FhoII = 6x10"^^ s-l andexlO-^'^ s"!. The long- (short-) dashed 
curves take FhcII to be 3 (10) times these values, but assume that 
equilibrium has not been reached and xhcII o = 1- The arrow- 
heads represent this Letter's constraint and an estimated future 
constraint, probably using Ly/3. The shaded region is A;, times the 
gas density PDF, arbitrarily normalized. 

From these two propositions, one can place an inter- 
esting lower bound on xhoII in ^ 10 cMpc regions. To see 
this, let us take our constraint of a;Heii > 0.008 at A;, — 
0.15 (Section [2]). Using the two propositions and equa- 
tion (4) to infer the ionization state of neighboring gas el- 
ements (which provides the scaling xucU oc 
we find that nearby elements with A^ = 1 should have 
a^Heii > 0.034. See the intersection of the solid verti- 
cal line in Figure [2] with the thick solid curve. In ad- 
dition, the intersection of this vertical line with the two 
dashed curves in Figure [2] represents the lower limit on 
the volume- and mass-weighted Hell fractions in the re- 
gion surrounding the opaque voids with ccHeii > 0.008 at 
Ab = 0.15, or xhcILV > 0.026 and xhoILM > 0.052 (ex- 
cluding Ah > 10 in these averages). This calculation of 
the volume- and mass-weighted Hell fractions uses the 
gas density probability distribution function (PDF) mea- 
sured from simulations in lMiralda-Escude et al.l ()2000[ ).^ 
and it uses the fact that a > 10 cMpc region is fairly 
representative of the IGM at z — 3.^ 

It is clear that these constraints can be improved with 
a more careful analysis and better data. We expect that 
the constraint a;Heii,v > 0.1 in a > 10 cMpc region is 
possible by utilizing smaller Af, (probably in a statistical 
manner) or by using higher Lyman-series resonances. 
For example, the measurement of r > 4 at A^ — 0.15 in 

* This PDF agrees well at 2 3 with more recent simula- 
tions wi th an updated cosmology in the range 0.1 < A;, < 10 
IIBolton fc Becker 2009). ^ ^ 

^ In linear theory in the assumed cosmology, the standard de- 
viation of A;, averaged over a sphere of radius 5 (10) cMpc is 0.4 
(0.2) at 2 = 3. 
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Ly/3 absorption would constrain xhcII.v > 0.1 (dashed 
vertical line in Figure 2). 

We have shown that if two propositions hold, the Hell 
Lyman forest can place strong constraints on xhoII in 
a region. Next, we argue that the two propositions are 
generally satisfied. 

4.1. Proposition 1: xyicU oc 0.7(7 1) 

We define Hell reionization as being over in a region 
when a;Hcii ^ 1 and T[t) is much greater than afif., 
such that the region never significantly recombines. The 
scaling xhoII oc 0.7(7 1) ^-^^l apply a few equilibrium 
times after HeH reionization ends in such a region even 
if photoionization equilibrium does not to a precision of 
(^xhoII < exp(— Ai(r}) (Equation [31), where (. . .} rep- 
resents a time average from the time HeH reionization 
ended. At. Furthermore, once this scaling is achieved 
after HeH reionization, this scaling is maintained even if 
photoionization equilibrium is not. 

During or a few icq after HeH reionization this scaling 
may not be achieved, but it is still difficult to signif- 
icantly evade the constraint on xhoII that assumes this 
scaling. To see this, let us consider a region with uniform 
rHoii(i) that has an initial HeH fraction of xhcII.o ~ 1 
and is being reionized (xhcII is decreasing). In this case, 
assuming photoionization equilibrium and using under- 
dense regions to infer the largest allowed photoionization 
rate Fg^f^ could result in an overestimate of the HeH 
fraction if the incident FhcII > Tg^j'^. This effect is il- 
lustrated by the long- and short-dashed curves in Figure 
[21 which are not equilibrium curves but which represent 
3 and 10 times larger FhoII than the corresponding solid 
curve (which assumes equilibrium). However, it is un- 
likely that a region will have FhcII ^ ^BsIii because the 
timescale to reach the new equilibrium state is ~ Tj^^jj, 
which is typically much less than H{z)~^ . For example, 
the constraint xhoII = 0.1 at Ab = 1 and z — i would 
imply icq < r^rfV' « QmH{z)-\ 

Figure [4l illustrates this argument for a toy example 
where FhcII is taken as constant after HeH reionization 
begins. Plotted are contours of fixed HeH fraction in the 
rHcii versus Atj-oi plane, where Airei is the time since 
the beginning of HeH reionization in the region. The 
top thick solid curve represents our constraint xhcII = 
0.008 at Ab = 0.15. Only the colored area left of this 
curve is allowed. The dashed curves represent a;Hcii v = 
0.01, 0.02, 0.04, 0.08, and 0.16. This illustrates that in 
order for xueii.v to be significantly less than the limit we 
set assuming photoionization equilibrium (or XHeii,v = 
0.026; Section 4), the HeH needs to have been recently 
reionized. For example, to evade this limit such that 
a^Hcii,v — 0.015, the HeH in the region must have been 
reionized over the last 100 Myr (Az w 0.1). 

4.2. Proposition 2: FhcII fluctuates on > 10 cMpc 

This proposition holds because the mfp of a HeH- 
ionizing photon, which sets the correlation length for 
FhcII, should be > lOcMpc. The mfp for a Hell Lyman- 
limit photon to be absorbed in diffuse gas at z = 3 
is w 0.8 a;]^^jj y cMpc (assuming a homogeneous IGM), 
which is > lOcMpc for the value of xhcII.v that satu- 
rates the previous bounds we have quoted, XHeii.v < 0.1. 
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Fig. 4. — Contours of fixed xhcH in the ruoll versus Airci plane 
at z = 3, where At^^i is the time since the beginning of Hell 
reionization in the region assuming constant FHcn- The upper 
solid curve represents this Letter's constraint of xjjcn = 0.008 at 
A(, = 0.15, and the colored region left of this curve is still allowed. 
The lower solid curve is the constraint that a future observation 
could achieve of XHeU = 0.04 at A;, = 0.15. The thick dashed 
curves represent a;Hcn,v = 0.01, 0.02, 0.04, 0.08, and 0.16. 

The spectrally-averaged mfp, weighting by contribution 
to FhcII, is 3 and 7 times larger than the Lyman- limit 
mfp for photon spectral indexes of 2.5 (expected in the 
case of no absorption) and 1.5, respectively. In addition, 
estimates for the mfp of He ll-ionizing photons to be ab- 
sorbed in dense systems a la lHaardt &: Madaul (|1996[ ) are 
simil ar at the FhsTT that y ield a;Hcii,cq.v < 0.1 (Appendix 
A in lMcQuinn et alll2009l plus foUowup work). In agree- 
ment wit h this proposition, the Hell reionization simu- 
lations of iMcQuinn et all ()2009l ) find FHoii-fluctuations 
on > 10 cMpc scales (even in Hell regions; see their Fig- 
ure 5). 

The previous paragraph argued that FhcII fluctuates 
as > 10 cMpc when our bound on a;Hcii.v is saturated. 
However, xhoII,v could be larger than our lower bound 
immediately around the underdensity such that the fluc- 
tuation scale is < 10 cMpc. One could then design a 
situation in which the average xhcII in a 10 cMpc region 
around an underdensity is less than our bound. How- 
ever, the situation would have to be very contrived in 
order to signiflcantly evade our bound, especially since 
current spectra show saturated voids throughout their 
Gunn-Peterson troughs. 

In appa r ent contradiction to this proposition, 
iShull et all (|2004l ) detected large fluctuations in FhcII 
on ~ 1 cMpc scales in the spectra of HE2347-4342, with 
the least den s e pixe ls yielding small er FHnTT. H o wever , 
iBolton et al] (pOOSh replicated the IShull et all ([200I ) 
analysis on mock data with the same noise properties, 
and they showed that the observed fluctuations could 
be explained by Poisson fluctuations in the distribution 
of quasars and by a 30cMpc mfp. It appears the 
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re ason that a ^ IcMpc correlation length is inferred 
m IBolton et all ((2005h owes to noise in the Hell forest 
spectrum, which often results in a la rge overestimate for 
FhcII in pixels with r > 0.05. In fact. lFechner fc ReimersI 
(poo 71) demonstrated the presence of this bias and, when 
it was corrected for, estimated a fluctuation scale for 
rHcii of > 8 - 24cMpc from HE2347-4342. 

Related to Proposition 2 is the misconception that 
FHeii should correlate strongly with the density of an ab- 
sorber. The rays from a quasar will traverse many (essen- 
tially uncorrelated) absorption systems prior to traveling 
~ 10 cMpc to a typical absorber, making it implausible 
that the typical overdense element will be exposed to a 
larger FHeii from these sources than an underdense one. 
This argument ignores that overdense absorbers are more 
likely to exist closer to a quasar where there is more ra- 
diation. However, the correlation between the locations 
of ~ quasars (~ 1 per 30'^ c Mpc'^) and the densit y 
of a typical absorber is negligible (jMcQuinn et al.|[2009f ). 
Even if quasars were not the primary source for 

THeii [contrary to what the quasar luminosity function 
and the amplitude of FHsTT- fluctuations in the HcH for- 
est suggest (McQu inn et al. 2009; Bolton ct al. 2005)], 
it would be difficult to make Fueii correlate strongly 
with the density of absorbers: as the number density 
of the sources increases, the fluctuations in FhcII should 
decrease. 

5. CONCLUSIONS 

This Letter showed that current Hell Lya forest data 
imply that shcII.v > 0.03 in regions surrounding opaque 



voids, while xhcILV ^0.1 should be possible with fu- 
ture HeH Lyman forest data. These limits derive from 
assuming photoionization equilibrium and that the max- 
imum allowed FhoII in an underdense gas element is also 
the maximum in a surrounding > 10 cMpc region - the 
expected fluctuation scale for FhcII- They may be weak- 
ened somewhat, but only if the region in question is be- 
ing reionized. However, since HeH reionization is the 
epoch when a non-negligible Hell fraction exists in the 
low-density IGM, by definition it would be occurring if 
a;Hcii.v ^0.1. Therefore, the HeH forest may be able to 
address whether HeH reionization is occurring at z sa 3. 

Our constraint of xhcII.v > 0.03 is a significant 
improvement over previous efforts, which set the bound 
a^HpTT.v > 1 . 3 X 1 0~^ in the regions of highest opacity 
([Heap et al.ll2000f l. Our constraint on xhcILV in select 
regions is 100 times stronger than the tightest constraint 
on the volume-averaged HI fraction from the z > 6 HI 
Lya-7 forests of xm,v > 3 x 10"^ (|Fan et all 120061) . It 
can be improved by better observations, by including 
redshift-space effects, by measuring xhoII in pixels with 
Af) < 0.15, or by utilizing absorption information from 
higher Lyman-series resonances. 
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